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Molecular  Orbital  Based  Design  Guideline  for  Hypergolic 
Ionic  Liquids 

Debasis  Sengupta*[a]  and  Ghanshyam  L.  Vaghjiani[b] 


Abstract:  Currently,  monomethyl  hydrazine  is  the  most 
widely  used  hypergolic  rocket  fuel.  Due  to  its  high  toxic 
vapor,  there  is  a  thrust  towards  developing  low-toxic  hyper¬ 
golic  fuels.  Ultra-low  vapor  pressure  ionic  liquids  are  one 
such  potential  category  of  fuels.  However,  designing  ionic 
liquid  with  ignition  delay  comparable  to  monomethyl  hy¬ 
drazine  is  a  challenge,  because  fundamental  understanding 
of  the  hypergolic  nature  of  ionic  liquids  is  far  from  clear. 
This  work  used  the  computed  energy  gap  values  between 
the  highest  occupied  molecular  orbitals  (HOMO)  of  the 
anions  for  a  series  of  ionic  liquids  and  the  lowest  occupied 
molecular  orbital  (LUMO)  of  HN03,  and  variation  in  the 
computed  relative  heats  of  formation,  A Hf,  of  these  anions 


to  develop  correlations  to  predict  hypergol  activity  be¬ 
tween  an  ionic  liquid  fuel  and  nitric  acid  as  the  oxidizer. 
The  observed  trends  in  HOMO-LUMO  energy  gap  and  A Hf 
values  can  be  used  successfully  to  verify  not  only  hypergo- 
licity  of  known  systems  but  also  the  lack  of  this  phenomen¬ 
on  in  OH"  and  BF4“  based  ionic  liquids.  It  was  shown  that 
through  suitable  substitution  of  electron  withdrawing  or 
electron  donating  groups  in  the  anion,  the  energy  gap  and 
the  AHf  values  could  be  tailored  into  an  optimal  range  that 
would  have  a  high  probability  for  the  new  system  to  exhib¬ 
it  hypergolic  reactivity.  To  validate  our  method,  we  suggest 
herein  new  ionic  liquid  structures  for  synthesis  and  experi¬ 
mental  screening. 


Keywords:  Hypergolic  •  Molecular  orbital  •  Design  guideline  •  Heat  of  formation  •  Energetic  ionic  liquids 


Introduction 

A  hypergolic  fuel  will  spontaneously  ignite  when  mixed 
with  a  suitable  oxidizer.  These  types  of  fuels  are  routinely 
used  as  rocket  propellants  in  many  space  and  missile  appli¬ 
cations.  Common  oxidizers  that  may  be  used  are  nitrogen 
tetroxide  (NTO)  and  inhibited  red  fuming  nitric  acid 
(IRFNA).  Currently  used  hypergolic  fuels  include  hydrazine, 
as  well  as  its  derivatives;  monomethyl  hydrazine  (MMH) 
and  unsymmetrical  dimethyl  hydrazine  (UDMH).  These  fuels 
are  highly  toxic,  have  high  vapor  pressures,  and  are  known 
to  be  carcinogens.  For  these  reasons,  replacements  for  hy¬ 
drazines  are  sought.  Ionic  liquids  (ILs),  due  to  their  negligi¬ 
ble  vapor  pressure,  are  being  considered  as  viable  alter¬ 
nates.  Since  the  recent  discovery  of  dicyanamide  (DCA) 
based  hypergolic  energetic  ionic  liquids  (EILs)  [1],  several 
papers  have  reported  on  the  synthesis  of  new  hypergolic 
EILs,  which  have  the  potential  to  provide  improved  propul¬ 
sion  performance  [2-4]  compared  to  the  MMH/NTO  base¬ 
line  system.  An  important  challenge  is  to  design  hypergolic 
EILs  with  short  ignition  delay  (ID)  times  that  are  compara¬ 
ble  to  the  ID  time  for  MMH/NTO.  Under  current  practices, 
EILs  are  first  synthesized  without  any  prior  knowledge  of 
their  ability  of  being  hypergolic  with  a  suitable  oxidizer.  A 
trial-and-error  drop  test  is  followed  to  discern  whether  the 
fuel  will  actually  undergo  a  hypergolic  reaction.  Therefore, 
the  development  of  a  simple,  but  reliable  design  guideline 
for  a  priori  prediction  of  hypergolic  EILs  could  have  a  signifi¬ 


cant  impact  in  reducing  the  time  to  develop  a  viable  hyper¬ 
golic  EIL  propellant. 

Although  there  have  been  some  efforts  in  the  past  to  un¬ 
derstand  hypergolic  reaction  mechanisms  and  to  correlate 
ID  times  with  molecular  properties  [5-8],  the  utility  of  such 
schemes  has  been  limited  to  molecular  liquids,  since  the  in¬ 
creased  complexity  of  cations/anions  Columbic  interactions 
present  in  ILs  is  not  required  to  be  considered  for  molecu¬ 
lar  liquids.  Recently,  we  reported  relationships  between  ID 
times  and  a  set  of  quantum  chemically  computed  descrip¬ 
tors  using  linear  quantitative  structure  property  relationship 
(QSPR)  and  artificial  neural  network  (ANN)  approaches  [9]. 
The  observed  correlations  were  based  on  charge  partial 
charge  surface  area  (CPSA)  based  descriptors  calculated 
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from  semi-empirical  quantum  calculations.  Chambreau 
et  al.  [10]  reported  Fourier  transform  infrared  experiments 
and  ab  initio  calculations  to  understand  the  mechanism  of 
hypergolic  reaction  of  a  DCA  anion  based  EIL  and  nitric 
acid,  and  found  the  underlying  reactions  to  be  very  com¬ 
plex.  We  have  performed  reactive  molecular  dynamics  sim¬ 
ulations  on  an  EIL/nitric  acid  system  to  understand  the  fun¬ 
damental  chemical  reactions  that  lead  to  hypergolic  igni¬ 
tion  [11]  and  also  observed  the  reaction  chemistry  to  be 
complex.  Although  these  studies  provide  important  infor¬ 
mation  regarding  the  initiation  chemistry,  there  is  no  guid¬ 
ance  available  in  the  literature  for  a  priori  design  of  ILs  that 
has  a  potential  to  show  hypergolicity.  Previous  experimen¬ 
tal  observations  have  concluded  that  the  anion  plays  an 
important  role  in  determining  the  hypergolic  tendency  of 
EILs  [1,10,12].  Recently,  McRay  et  al.  [13]  attempted  to  cor¬ 
relate  ID  times  with  molecular  structures  and  thermophysi¬ 
cal  properties,  such  as  viscosity,  thermal  conductivity  and 
density.  While  they  found  some  evidence  that  ID  times  de¬ 
pended  on  the  nature  of  the  exocyclic  substitution  in  the 
five-membered  ring  cations  for  dicyanamide-based  ILs,  no 
significant  correlation  was  observed  with  any  of  the  ther¬ 
mophysical  properties. 

The  goal  of  this  work  is  to  obtain  a  fundamental  under¬ 
standing  of  the  cause  of  hypergolic  reactions  based  on  the¬ 
oretical  principles  and  provide  a  guideline  for  designing  hy¬ 
pergolic  ionic  liquids.  In  principle,  such  a  guideline  may 
therefore  be  used  by  a  synthetic  chemist  at  the  concept 
stage  for  selecting  an  anion,  which  whould  have  a  high 
probability  for  forming  a  hypergolic  fuel.  The  validity  of  the 
present  approach  is  tested  against  a  number  of  anions  that 
have  been  reported  in  the  literature. 


Computational  Methods 

Quantum  mechanical  calculations  were  performed  using 
B3LYP  exchange  correlation  functional  [14-16]  with  6-31  + 
+  G(d,p)  basis  set.  Diffuse  functions  were  added  to  obtain 
a  better  description  of  the  anions.  All  heats  of  formation 
values,  AHf,  reported  here  are  calculated  at  298.15  K  using 
a  highly  accurate  complete  basis  set  method,  CBS-QB3  [17]. 
The  mean  absolute  deviation  of  the  heats  of  formation  of 
the  method  is  reported  to  be  4.98  Idmol-1  [18].  All  calcula¬ 
tions  were  performed  using  Gaussian  03  program  package 
[19]. 

Results  and  Discussion 

Importance  of  the  Anion 

Figure  1  shows  the  anions  that  were  experimentally  synthe¬ 
sized  and  the  hypergolicity  of  the  EILs  containing  these 
was  experimentally  tested. 

Figure  2  shows  that  the  nature  of  the  electron  density 
surface  of  the  HOMO  of  the  cation/anion  ion-pair  remains 
unchanged  compared  to  that  of  the  isolated  anion,  and  ac- 
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Figure  1.  Hypergolic  EILs,  whose  ignition  delay  (ID)  times  were  re¬ 
ported  (see  for  example  Refs.  [1,11]). 


Figure  2.  HOMO  of  the  ion-pair  (upper  panel)  compared  to  that  of  the  corresponding  isolated  anion  (lower  panel). 
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cording  to  frontier  molecular  orbital  (FMO)  theory  [20],  the 
HOMO  of  the  EIL  should  be  an  electron  donor,  where  the 
reaction  between  the  EIL  and  the  nitric  acid  can  be  consid¬ 
ered  to  be  of  the  acid-base  type.  Previously,  this  FMO  ap¬ 
proach  has  been  applied  in  interpreting  reactivity  in  Diels— 
Alder  reactions,  photochemical  reactions,  and  organometal- 
lic  reactions  [21-23]. 

According  to  the  molecular  orbital  representation  of  the 
acid-base  theory,  the  interaction  between  the  acid  and  the 
base  is  dictated  by  the  "energy  gap"  between  the  HOMO 
of  the  anion  and  LUMO  of  HN03,  and  this  quantity  can  be 
explored  as  a  criterion  for  hypergolicity  (or  reactivity).  While 
it  is  mathematically  complex  to  construct  an  accurate  inter¬ 
action  matrix  element  of  the  orbital  interaction,  it  is  known 
from  perturbation  theory  that  the  orbital  interaction  is  dic¬ 
tated  by  the  energy  gap  between  orbitals  [24];  the  smaller 
the  gap,  the  larger  the  interaction.  Figure  3  shows  the 
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Figure  3.  Computed  "energy  gap"  for  various  experimentally  stud¬ 
ied  anions  (see  Tables  SI  and  S2,  Supporting  Information,  for  the 
corresponding  heat  of  formation  and  ionization  potential  values). 

energy  gap  values  for  the  anions  studied  herein  and  for 
which  experimental  hypergolicity  is  known.  DCA  based  hy- 
pergolic  EILs  have  received  the  greatest  attention  in  the  lit¬ 
erature  following  the  first  discovery  of  an  IL  hypergol 
[1,12].  Therefore,  here  we  consider  DCA  as  the  baseline 
system  (horizontal  line  in  Figure  3).  Any  anion  with  an 
energy  gap  that  falls  in  the  vicinity  or  below  this  line  is  pre¬ 
dicted  to  be  hypergolic  when  coupled  with  an  appropriate 
cation  provided  it  meets  the  heat  of  formation  criterion  (as 
discussed  later). 

One  of  the  most  striking  features  is  the  difference  in  the 
"energy  gap"  values  of  azide  and  DCB.  Azide  ILs  reported  in 
the  literature  were  found  to  be  non-hypergolic,  except  for 
two  examples  [12].  On  the  other  hand,  DCB  containing  EILs 
have  been  found  to  be  hypergolic.  EILs  containing  NCA 
were  also  found  to  be  hypergolic.  N03_  has  a  higher 
"energy  gap"  than  DCA  and  is  predicted  to  be  non-hyper- 


golic.  In  fact,  so  far  only  one  IL  containing  N03~  was  found 
to  be  hypergolic.  Recently,  a  few  EILs  containing  the 
H2P02_  (HP)  anion  have  been  synthesized  by  the  Shreeve 
group  [25].  The  majority  of  them  were  found  to  be  hyper¬ 
golic,  however,  they  had  longer  ID  times.  Also,  out  of  two 
synthesized  DCF  containing  EILs,  one  of  them  was  found  to 
be  hypergolic.  While  there  are  many  reports  of  EILs  con¬ 
taining  NCA,  DCA,  DCB,  and  azide  in  the  literature,  more 
H2P02-  and  DCF  containing  EILs  should  be  synthesized  to 
further  validate  our  probabilistic  approach  for  identifying 
EIL  hypergols. 

DMP  is  one  anion  that  has  a  lower  "energy  gap",  but  the 
only  DMP  EIL  so  far  synthesized  has  been  found  to  be  non- 
hypergolic.  Similarly,  two  other  anions,  vinyl  analogue  of 
DCA  (i.e.,  vDCA)  and  CF  [cyanoform,  C(CN)3_]  are  also  pre¬ 
dicted  to  be  hypergolic,  but  the  synthesized  EILs  contain¬ 
ing  these  two  anions  were  also  found  to  be  non-hypergol¬ 
ic.  A  close  look  at  the  structures  of  the  anions  reveals  that, 
unlike  other  anions,  both  vDCA  and  CF  contain  C— C  bonds. 
The  structural  difference  of  these  two  anions  from  others 
could  lead  to  differences  in  the  reactivity  that  is  not  cap¬ 
tured  in  the  simple  interpretation  described  above.  There¬ 
fore,  more  theoretical  and  experimental  work  needs  to  be 
performed  to  isolate  the  effect  of  C-C  bond(s)  over  C-N 
bond(s).  Although  it  has  been  shown  that  the  nature  of  the 
electron  density  surface  of  the  HOMO  of  the  cation/anion 
ion-pair  remains  unchanged  compared  to  that  of  the  isolat¬ 
ed  anion,  it  is  worth  discussing  the  effect  of  the  cation  on 
the  "energy  gap".  In  the  ion-pair  system,  there  is  no  cova¬ 
lent  bond  and  therefore  orbital  mixing  is  negligible.  How¬ 
ever,  orbital  energies  of  anions  are  significantly  affected  by 
the  electrostatic  field  created  by  the  positively  charged 
cation  (and  vice  versa).  This  electrostatic  effect  results  in 
stabilization  of  the  anion  HOMO.  Our  calculations  show 
that  for  a  variety  of  cations,  the  magnitude  of  this  stabiliza¬ 
tion  is  near  identical  for  a  particular  anion  and  that  the 
HOMO  energy  levels  of  the  ion-pairs  are  virtually  independ¬ 
ent  of  the  nature  of  the  cations.  Table  S4  (Supporting  Infor¬ 
mation)  shows  the  variance  in  HOMO  energy  level  is  below 
0.04  eV  for  the  azide  anion.  Since  the  LUMO  of  the  HN03  is 
constant,  one  could  calibrate  the  "energy  gap"  to  arrive  at 
a  new  rule.  However,  obtaining  optimized  ion-pair  struc¬ 
tures  is  a  non-trivial  task  due  to  the  non-bonded  interac¬ 
tions,  and  therefore  this  work  has  primarily  focused  just  on 
the  HOMO  of  the  anions. 

Anion  Design  Guideline 

Based  on  the  theoretical  ideas  discussed  above,  the  anion 
appears  to  play  an  important  role  in  EIL  hypergolicity.  The 
interpretation  is  based  on  the  energies  of  the  HOMO  of 
anions  relative  to  the  LUMO  of  the  oxidizer.  In  principle, 
one  can  engineer  the  HOMO  energy  of  an  anion  by  substi¬ 
tution  of  a  parent  anion  with  different  functional  groups.  It 
is  known  that  electron  withdrawing  groups  lower  the 
HOMO  energy,  while  electron  donating  groups  elevate  the 
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Figure  4.  Values  of  the  energy  gap  for  some  conceptual  anions. 
DCA  and  DCB  are  also  shown  for  comparison. 


energy.  Therefore,  — F,  — N02,  and  — CN  will  lower,  and  alkyl, 
amino,  and  substituted  amino  will  elevate  the  HOMO 
energy  compared  to  its  unsubstituted  analogue.  Since  a  hy¬ 
pergolic  reaction  ultimately  releases  a  large  amount  of 
heat,  the  energy  content  of  the  anion  must  also  be  taken 
into  account  in  the  design  strategy;  substitution  should  be 
curried  out  in  such  u  way  that  does  not  reduce  the  energy 
content  (i.e.,  the  A Hf)  of  the  new  anion  drastically  compared 
to  that  of  the  parent  anion.  Figure  4  shows  some  conceptu¬ 
al  anions  and  the  energy  gaps  between  their  HOMOs  and 
the  LUMO  of  the  HN03.  Their  A Hf  values  are  depicted  in 
Figure  5. 

Important  observations  to  note  (see  Table  S3,  Supporting 
Information)  include  the  following:  substituting  one  H  in 
DCB  by  -N02,  reduces  its  A Hf  from  -61.5  kJ  mol-1  to 
-214.4  kJ  mol-1,  while  a  -F  substituted  DCB  has  A Hf  of 
-418,68  kJ  mol-1.  Also,  -N02  brings  the  HOMO  energy  of 
the  anion  closer  to  the  LUMO  than  does  — F.  Therefore,  — 
N02  may  be  preferable  to  — F,  if  a  lowering  of  anion  HOMO 
is  required.  If  all  four  groups  of  DCB  are  replaced  by  -F,  the 
HOMO  of  the  resulting  BF4-  becomes  much  lower  than 
LUMO  of  HN03  resulting  in  a  lager  absolute  gap  than  with 
DCA  (i.e.,  1.54  eV).  Moreover,  its  A Hf  is  drastically  reduced 
(i.e.,  -1728.3  kJ  mol-1).  Therefore  the  present  theory  pre¬ 
dicts  BF4-  is  unlikely  to  be  hypergolic.  Indeed,  BF4-  contain¬ 
ing  ionic  liquids  are  found  to  be  very  stable,  non-energetic, 
and  non-hypergolic.  Therefore,  a  trade-off  between  the 
"energy  gap"  value  and  the  change  in  A Hf  should  be  con¬ 
sidered  during  any  substitution  strategy.  So  far  in  the  litera¬ 
ture,  HP  (with  energy  gap  similar  to  DCA)  is  the  lowest  A Hf 
anion  (—661.9  kJmol-1)  for  which  EIL  hypergolicity  was  ob¬ 
served.  Using  this  A Hf  anion  value  as  another  baseline  limit, 
we  see  that  the  A Hf  of  DMP  of  -1128.4  kJ  mol-1  is  far 
below  this  limit  and  has  a  greater  influece  than  its  low 
"energy  gap"  of  Figure  2,  which  makes  this  IL  non-hyper- 
golic.  New  EILs  containing  DMP  should  be  synthesized  and 
tested  for  hypergolicity. 
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Figure  5.  A Hf  values  of  some  conceptual  anions.  DCA  and  HP,  re¬ 
spectively,  have  highest  and  lowest  values  of  anion  A Hf,  for  which 
EIL  hypergolicity  is  experimentally  observed. 


Promising  conceptual  anions  could  be  formed  via  a  sub¬ 
stitution  of  one  of  the  H's  of  DCB  by  -N3.  Single  substitu¬ 
tion  would  slightly  elevate  the  HOMO  of  the  DCB,  but  the 
computed  A  Hf  increases  drastically,  from  -61,5  to 
158.7  kJmol-1.  Substitution  of  the  second  H  by  the  electron 
donating  -N(H)NH(CH3),  which  is  derived  from  MMH, 
would  elevate  the  DCB  HOMO,  but  would  still  remain 
below  the  baseline  value  of  DCA.  The  computed  A Hf,  how¬ 
ever,  would  increase  to  233.6  kJ  mol-1.  The  other  example  is 
the  substitution  of  one  -CN  group  of  DCA  by  -N02.  EILs 
with  this  anion  (NCA)  have  been  studied  and  in  many  cases 
found  to  be  hypergolic  [3].  This  substitution  brings  the 
HOMO  closer  to  HN03  LUMO  with  an  acceptable  reduction 
of  A Hf  (from  118.5  kJmol-1  for  DCA  to  -27.2  kJmol-1  for 
NCA).  For  the  OH-  anion,  our  calculations  show  that  the 
value  of  the  "energy  gap"  is  4.54  eV,  which  is  much  higher 
than  that  of  the  azide  anion,  and  with  A Hf  of 
-130.2  kJ  mol-1.  Clearly  OH-  anion  based  EILs  are  unlikely 
to  be  hypergolic. 

Finally,  based  on  the  above  computational  results  and 
validations,  we  provide  simple  rules  for  designing  anions 
(see  Figure  6  above):  (1)  Conceptualize  a  preliminary  struc¬ 
ture  for  an  anion  avoiding  any  C-C  bonds;  (2)  compute  its 
A Hf,  (3)  calculate  the  "energy  gap",  i.e.,  the  difference  in  en¬ 
ergies  of  anion  HOMO  and  LUMO  of  HN03;  (4)  compare 
with  the  baseline  value  (i.e.  of  DCA)  and  adjust  this  energy 
difference  by  functional  group  substitution,  keeping  the 
anion  A Hf  in  the  vicinity  of  other  hypergolic  anions  known 
experimentally;  and  (5)  combine  the  anion  with  a  high 
energy  content  cation  to  make  a  hypergolic  EIL.  It  should 
be  noted  that  the  theory  presented  here  does  not  consider 
the  effect  of  thermophysical  properties  and  therefore,  does 
not  predict  properties  such  as  melting  points,  diffusion  co¬ 
efficients,  viscosity,  thermal  conductivities,  and  specific 
heats  of  the  suggested  ILs  containing  the  conceptual 
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anions.  These  thermophysical  properties  will  control  the 
mixing  and  heat  transfer,  and  will  ultimately  affect  the  ID 
times  for  ILs  with  similar  chemical  reactivities.  In  order  to 
take  into  account  these  effects,  one  needs  to  first  compute 
these  properties  using  other  theoretical  methods,  such  as 
molecular  dynamics  simulations. 


Conclusions 

This  work  provides  a  fundamental  understanding  of  hyper- 
golic  reaction  between  energetic  ionic  liquids  and  nitric 
acid  by  coupling  frontier  molecular  orbital  theory  and  ther¬ 
mochemistry  derived  from  ab  initio  quantum  chemistry.  A 
simple  theoretical  model  based  on  frontier  molecular  orbi¬ 
tal  theory  was  proposed  for  designing  anions.  The  HOMO 
of  the  anion  retains  its  character  in  the  ion-pair  and  there¬ 
fore  expected  to  be  the  dominant  reactant.  The  energy 
gap  between  the  HOMO  of  anions  and  the  LUMO  of  HNOs 
was  used  as  a  parameter  for  interaction  between  a  base 
and  an  acid.  The  larger  the  interaction  is,  the  higher  is  the 
probability  that  the  EIL  will  be  hypergolic  provided  it  meets 
the  requirement  of  heat  of  formation.  Several  concept 
anions  are  proposed  here  that  have  the  potential  to  form 
hypergolic  EILs.  Several  other  factors,  such  as  mixing  rate  of 
the  EIL  and  the  oxidizer,  heat  transfer  rates,  reaction  kinet¬ 
ics  and  the  chemical  influence  of  the  cation  are  capable  of 
affecting  hypergolicity,  and  the  scope  of  the  present  design 
approach  could  be  expanded  in  the  future  to  fully  explore 
these  effects. 

Supporting  Information  (see  footnote  on  the  first  page 
of  this  article):  Details  of  the  heat  of  formation  and  energy 
gap  values. 
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